-- declaring species
species = {
    { "H+",     1 },
    { "HO-",   -1 },
    { "OxH2",   0 },
    { "OxH-",  -1},
    { "Ox--",  -2 },
    { "Na+",    1 },
    { "Cl-",   -1 },
    { "NH3",    0 },
    { "NH4+",   1 }
};

eqs =
{
    { "water",    10^(-14),  {1,"H+"},  {1,"HO-"} },
    { "Ox1",      10^(-1.27), {1,"OxH-"},  {1,"H+"}, {-1, "OxH2"}  },
    { "Ox2",      10^(-4.27), {1,"Ox--"},  {1,"H+"}, {-1, "OxH-"}  },
    { "ammoniac", 10^(-9.2), {1,"NH3"}, {1,"H+"}, {-1, "NH4+"}}
};

C0 = 0.01;
ini =
{
    { "E/N" },
    { 0, {1,"Na+"} },
    { 0, {1,"Cl-"} },
    { C0, {1,"OxH2"}, {1,"OxH-","Ox--"} },
    { 0,  {1,"NH3"}, {1,"NH4+"} }
};

C1 = 0.01;
dose =
{
    { "E/N" },
    { C0, {1,"Na+"} },
    { 0,  {1,"Cl-"} },
    { 0,  {1,"OxH2"}, {1,"OxH-","Ox--"} },
    { 0,  {1,"NH3"}, {1,"NH4+"} }
}

V0  = 10;
V1  = V0*3;